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Abstract 

We compute the big bang nucleosynthesis limit on the number of light neutrino degrees 
of freedom in a model-independent likelihood analysis based on the abundances of 4 He 
and 7 Li. We use the two-dimensional likelihood functions to simultaneously constrain 
the baryon-to-photon ratio and the number of light neutrinos for a range of 4 He 
abundances Y p = 0.225 - 0.250, as well as a range in primordial 7 Li abundances from 
(1.6 to 4.1) xl0~ 10 . For ( 7 Li/H) p = 1.6 x 10~ 10 , as can be inferred from the 7 Li 
data from Population II halo stars, the upper limit to N v based on the current best 
estimate of the primordial 4 He abundance of Y p = 0.238, is N v < 4.3 and varies from 
N v < 3.3 (at 95% C.L.) when Y p = 0.225 to N v < 5.3 when Y v = 0.250. If 7 Li is 
depleted in these stars the upper limit to N v is relaxed. Taking ( 7 Li/H) p = 4.1 x 10~ 10 , 
the limit varies from N v < 3.9 when Y p = 0.225 to N v < 6 when Y p = 0.250. We also 
consider the consequences on the upper limit to N v if recent observations of deuterium 
in high-redshift quasar absorption-line systems are confirmed. 
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1 Introduction 



One of the most important limits on particle properties is the limit on the number of light 
particle degrees of freedom at the the time of big bang nucleosynthesis (BBN) [|T| . This is 
commonly computed as a limit on the number of light neutrino flavors, N u . Recently, we 
used a model-independent likelihood method (see also || |[]) to simultaneously constrain 
the value of the one true parameter in standard BBN, the baryon-to-photon ratio 77, together 
with N u . For similar approaches, see |J. In that work ||, we based our results on the best 
estimate of the observationally determined abundance of 4 He, Y p = 0.234 ± 0.002 ± 0.005 
from ||, and of 7 Li, 7 Li/H = (1.6 ± 0.1) x 10 -10 , from [0. While these determinations can 
still be considered good ones today, there is often discussion of higher abundance for 4 He as 
perhaps indicated by the data of || and higher abundances of 7 Li due to the effects of stellar 
depletion (see e.g. ||). Rather than be forced to continually update the limit on N v as the 
observational situation evolves, we generalize our previous work here and compute the upper 
limit on N v for a wide range of possible observed abundances of 4 He and 7 Li. Because the 
determinations of D/H in quasar absorption system has not dramatically improved, we can 
only comment on the implications of either the high or low D/H measurements. 

One of the major obstacles in testing BBN using the observed abundances of the light 
element isotopes rests on our ability to infer from these observations a primordial abundance. 
Because 4 He, in extragalactic HII regions, and 7 Li, in the atmospheres of old halo dwarf stars, 
are both measured in very low metallicity systems (down to l/50th solar for 4 He and l/1000th 
solar for 7 Li), very little modeling in the way of galactic chemical evolution is required to 
extract a primordial abundance for these isotopes. Of course systematic uncertainties, such 
as underlying stellar absorption, in determining the 4 He abundance and the effects of stellar 
depletion of 7 Li lead to uncertainties in the primordial abundances of these isotopes, and it is 
for that reason we are re-examining the limits to N v . Nevertheless, the problems in extracting 
a primordial 4 He and 7 Li abundance pale in comparison with those for D and 3 He, both of 
which are subject to considerable uncertainties not only tied to the observations, but to 
galactic chemical evolution. In fact, 3 He also suffers from serious uncertainties concerning 
its fate in low mass stars |10j. 3 He is both produced and destroyed in stars making the 
connection to BBN very difficult. 

Deuterium is totally destroyed in the star formation process. As such, the present or 
solar abundance of D/H is highly dependent on the details of a chemical evolution model, 
and in particular the galactic star formation rate. Unfortunately, it is very difficult at the 
present time to gain insight on the primordial abundance of D/H from chemical evolution 
given present and solar abundances since reasonably successful models of chemical evolution 
can be constructed for primordial D/H values which differ by nearly an order of magnitude^ 
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Of course much of the recent excitement surrounding deuterium concerns the observation 
of D/H in quasar absorption systems [n]-[r7|. If a single value for the D/H abundance in 

1 There may be some indication from studies of the luminosity density at high redshift which implies 
a steeply decreasing star formation rate [Q, and that at least on a cosmic scale, significant amounts of 
deuterium has been destroyed O]. 
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these systems could be established^], then one could avoid all of the complications concerning 
D/H and chemical evolution, and because of the steep monotonic dependence of D/H on r], a 
good measurement of D/H would alone be sufficient to determine the value of 77 (since D/H 
is nearly independent of N v ). In this case, the data from 4 He and 7 Li would be most valuable 
as a consistency test on BBN and in the case of He, to set limits on particle properties. In 
the analysis that follows, we will discuss the consequences of the validity of either the high 
or low D/H determinations. 

Using a likelihood analysis based on 4 He and 7 Li ||], a probable range for the baryon-to- 
photon ratio, 77 was determined. The 4 He likelihood distribution has a single peak due to the 
monotonic dependence of 4 He on 77. However, because the dependence on 77 is relatively flat, 
particularly at higher values of Y p , this peak may be very broad, yielding little information 
on 77 alone. On the other hand, because 7 Li is not monotonic in 77, the BBN prediction has 
a minimum at r] w ~ 3 (r7 10 = 10 10 r]) and as a result, for an observationally determined value 
of 7 Li above the minimum, the 7 Li likelihood distribution will show two peaks. The total 
likelihood distribution based on 4 He and 7 Li is simply the product of the two individual 
distributions. In M, the best fit value for 7710 based on the quoted observational abundances 
was found to be 1.8 with a 95% CL range 

1.4 < Vl0 < 4.3 (1) 

when restricting the analysis to the standard model, including N u = 3. In determining (|TJ) 
systematic errors were treated as Gaussian distributed. When D/H from quasar absorption 



systems (those showing a high value for D/H [14, jlBj ) is included in the analysis this range 
is cut to 1.50 < 7710 < 2.55. 

In ||, the maximum likelihood analysis of ||, |J which utilized a likelihood function C(rj) 
for fixed N v = 3 was generalized to allow for variability in N v . There a more general likelihood 
function £(77, N u ) was applied to the current best estimates of the primordial 4 He and 7 Li 
abundances. Based on the analysis in M, we chose Y p = 0.234 ± 0.002(stat.) ± 0.005(syst.) 
as well as the lower value Y p = 0.230 ± 0.003(stat.) ± 0.005(syst.) based on a low metallicity 
subset of the data. Using these values of Y p along with the value (Li/H) p = (1.6±0.07) x 10~ 10 
from [7J], we found peak likelihood values 7? 10 = 1.8 and N v = 3.0 with a 95% CL range of 
1.6 < N v < 4.0, 1.3 < 7710 < 5.0 for the higher 4 He value and similar results for the lower 
one. More recent data from Izotov and Thuan |19j seems to indicate a still higher value for 



Y p , and for this reason as well as wishing to be independent of the "current' best estimate of 
the abundances, we derive our results for a wide range of possible values for Y p and (Li/H) p 
which will account for the possibility of stellar depletion for the latter ||. Finally, in [[|, 
we considered only the effect of the high D/H value from quasar absorption systems. Since 
there was virtually no overlap between the likelihood functions based on the low D/H value 
and the other two elements, there was little point in using that value in our analysis. Since 
then, the low D/H value has been raised somewhat, and that together with our present 
consideration of higher Y p and (Li/H) p values makes the exercise worth while. 

2 It is not possible that all disparate determinations of D/H represent an inhomogeneous primordial 
abundance as the corresponding inhomogeneity in r\ would lead to anisotropics in the microwave background 
in excess of those observed fn8| . 
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In this paper, we follow the approach of M in constraining the theory on the basis 
of the 4 He and 7 Li data and to a lesser extent D/H, by constructing a likelihood function 
J~-(rj, N u ). We discuss the current status of the data in section 2, and indicate what range 
of values for the primordial abundances we consider. In section 3, we display the likelihood 
functions we use. As this was discussed in more detail in 0, §], we will be brief here. Our 
results are given in section 4, and we draw conclusions in section 5. 



2 Observational Data 

Data pertinent to the primordial 4 He abundance is obtained from observations of extragalac- 
tic HII regions. These regions have low metallicities (as low as l/50th solar), and thus are 
presumably more primitive than similar regions in our own Galaxy. The 4 He abundance 
used to extract a primordial value spans roughly an order of magnitude in metallicity (e.g. 
O/H). Furthermore, since there have been a considerable number of such systems observed 
with metallicities significantly below solar, modeling plays a relatively unimportant role in 
obtaining the primordial abundance of 4 He (see e.g. ||20|| ). 



The 4 He data based on observations in |2T|, § were discussed in detail in ||. There 
are over 70 such regions observed with metallicities ranging from about 2-30% of solar 
metallicity. This data led to the determination of a primordial 4 He abundance of Y p = 
0.234 ± 0.002(stat.) ± 0.005(syst.) used in 0. That the statistical error is small is due to the 
large number of regions observed and to the fact that the 4 He abundance in these regions 
is found to be very well correlated to metallicity. In fact, as can be understood from the 
remarks which follow, the primordial 4 He abundance is dominated by systematic rather than 
statistical uncertainties. 

The compilation in || included the data of ||. Although this data is found to be 
consistent with other data on a point by point basis, taken alone, it would imply a somewhat 
higher primordial 4 He abundance. Furthermore, the resulting value of Y p depends on the 
method of data analysis. When only 4 He data is used to self-consistently determine the 4 He 
abundance (as opposed to using other data such as oxygen and sulphur to determine the 
parameters which characterize the HII region and are needed to convert an observation of a 
4 He line strength into an abundance), a value of Y p as high as 0.244 ± 0.002 can be foundP] 

The problem concerning 4 He has been accentuated recently with new data from Izotov 



and Thuan [19]. The enlarged data set from [21,119] was considered in [2C[. The new resulting 



value for Y p is 



Y p = 0.238 ± 0.002(stat.) ± 0.005(syst.) (2) 



The new data taken alone gives Y p = 0.2444 ± 0.0015 when using the method based on a 
set of 5 helium recombination lines to determine all of the H II region parameters. By more 
conventional methods, the same data gives Y p = 0.239 ±0.002. As one can see, the 4 He data 
is clearly dominated by systematic uncertainties. 



3 We note that this method has been criticized as it relies on some 4 He data which is particularly uncertain, 
and these uncertainties have not been carried over into the error budget in the 4 He abundance fl6|. 
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There has been considerably less variability in the 7 Li data over the last several years. 
The 7 Li abundance is determined by the observation of Li in the atmospheres of old halo 
dwarf stars as a function of metallicity (in practice, the Fe abundance). The abundance used 
in [[| from the work in |7j continues to lead to the best estimate of the 7 Li abundance in the 
so called Spite plateau 

y 7 = ^ = (1.6 ± 0.07) xlO- 10 (3) 

where the error is statistical, again due to the large number of stars observed. If we employ 
the basic chemical evolution conclusion that metals increase linearly with time, we may infer 
this value to be indicative of the primordial Li abundance. 

In 0, we noted that there are considerable systematic uncertainties in the plateau abun- 
dance. It is often questioned as to whether the Pop II stars have preserved their initial 
abundance of Li. While the detection of the more fragile isotope 6 Li in two of these stars 
may argue against a strong depletion of 7 Li p2 , Q , it is difficult to exclude depletion of the 
order of a factor of two. Therefore it seems appropriate to allow for a wider range in 7 Li 
abundances in our likelihood analysis than was done in ||. 

There has been some, albeit small, change in the D/H data from quasar absorption 
systems. Although the re-observation of the high D/H in |23] has been withdrawn, the 
original measurements (14| of this object still stand at the high value. More recently, a 
different system at the relatively low redshift of z = 0.7 was observed to yield a similar high 
value [TTJ 

y 2 = D/H = (2 ±0.5) x 10" 4 . (4) 



The low values of D/H in other such systems reported in [15] have since been refined to show 
slightly higher D/H values [17] 

y 2 = D/H = (3.4 ±0.3) x 10" 5 . (5) 

Though this value is still significantly lower than the high D/H value quoted above, the low 
value is now high enough that it contains sufficient overlap with the ranges of the other light 
elements considered to warrant its inclusion in our analysis. 



3 Likelihood Functions 

Monte Carlo and likelihood analyses have been discussed at great length in the context of 
BBN g§ |25], [21| 0, [2§ §, ]§ 0. Since our likelihood analysis follows that of § and [g|, 
we will be very brief here. The likelihood function for 4 He, L 4 (iVy,?7) is determined from a 
convolution of a theory function 

(where Y p (N u ,r]) and aY{N u ,i]) represent the results of the theoretical calculation) and an 
observational function 

L.,oUY) = -J— exp f "^"^ 2 ) (7) 
V27rcr yo V 2(J Y0 J 
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where Yq and ayo characterize the observed distribution and are taken from Eqs. (0) and 
The full likelihood function for 4 He is then given by 

L,(N U , rj)= JdY L 4)Theory (F, N v , r,)L At0hs {Y) (8) 
which can be integrated (assuming Gaussian errors as we have done) to give 



(9) 



The likelihood functions for 7 Li and D are constructed in a similar manner. The quantities 
of interest in constraining the N v — r\ plane are the combined likelihood functions 

£47 = L 4 x L 7 (10) 

and 

L247 = L>2 x L47. (11) 

Contours of constant L 47 (or L247 when we include D in the analysis) represent equally likely 
points in the N v -rj plane. Calculating the contour containing 95% of the volume under the 
L 47 surface gives us the 95% likelihood region. From these contours we can then read off 
ranges of N u and rj. 



4 Results 

Using the abundances in eqs (§,^|) and adding the systematic errors to the statistical errors 
in quadrature we have a maximum likelihood distribution, L 47 , which is shown in Figure la. 
This is very similar to our previous result based on the slightly lower value of Y p . As one can 
see, L47 is double peaked. This is due to the minimum in the predicted lithium abundance 
as a function of 77, as was discussed earlier. We also show in Figures lb and lc, the resulting 
likelihood distributions, L 247 , when the high and low D/H values from Eqs. (^) and (|5p are 
included. 

The peaks of the distribution as well as the allowed ranges of rj and N v are more easily 
discerned in the contour plots of Figure 2 which shows the 50%, 68% and 95% confidence 
level contours in L47 and L247. The crosses show the location of the peaks of the likelihood 
functions. Note that L 47 peaks at N v = 3.2, (up slightly from the case with Y p = .234) and 
r] 10 = 1.85. The second peak of L 47 occurs at N u = 2.6, r] w = 3.6. The 95% confidence level 
allows the following ranges in rj and N v 

1.7 < N v < 4.3 

1-4 < r?i < 4.9 (12) 
These results differ only slight from those in 0. 
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Since L 2 picks out a small range of values of 77, largely independent of N u , its effect on 
L247 is to eliminate one of the two peaks in L 47 . With the high D/H value, L 2 47 peaks at 
the slightly higher value N v = 3.3, 7710 = 1.85. In this case the 95% contour gives the ranges 

2.2 <N U < 4.4 

1.4 < 7710 < 2.4 (13) 

(Strictly speaking, r] w can also be in the range 3.2 — 3.5, with 2.5 < N v < 2.9 as can be 
seen by the 95% contour in Figure 2a. However this "peak" is almost completely invisible 
in Figure lb.) The 95% CL ranges in N v for both L 47 and L 2 47 include values below the 
canonical value N v = 3. Since one could argue that N v > 3, we could use this condition as 
a Bayesian prior. This was done in |3C] and in the present context in B. In the latter, the 



effect on the limit to N„ was minor, and we do not repeat this analysis here. 

In the case of low D/H, L 2 picks out a smaller value of N v = 2.4 and a larger value of 
77 = 4.55. The 95% CL upper limit is now N v < 3.2, and the range for 77 is 3.9 < 7710 < 5.4. 



It is important to stress that with the increase in the determined value of D /H [17| in the 
low D/H systems, these abundances are now consistent with the standard model value of 
N v = 3 at the 2 a level. 

Although we feel that the above set of values represents the current best choices for the 
observational parameters, our real goal in this paper is to generalize these results for a wide 
range of possible primordial abundances. To begin with, we will fix (Li/H) p from Eq. (|3|), 
and allow Y v to vary from 0.225 - 0.250. In Figure 3, the positions of the two peaks of 
the likelihood function, L47, are shown as functions of Y p . The low-77 peak is shown by the 
dashed curve, while the high-77 peak is shown as dotted. The preferred value of N v = 3, 
corresponds to a peak of the likelihood function either at Y p = 0.234 at low 7710 = 1.8 or 
at Yp = 0.243 at 7710 = 3.6 (very close to the value of Y p quoted in |19[]). Since the peaks 
of the likelihood function are of comparable height, no useful statistical information can be 
extracted concerning the relative likelihood of the two peaks. The 95% CL upper limit to 
N v as a function of Y p is shown by the solid curve, and over the range in Y p considered varies 
from 3.3 - 5.3. The fact that the peak value of N v (and its upper limit) increases with Y p 
is easy to understand. The BBN production of 4 He increases with increasing N v . Thus for 
fixed Li/H, or fixed 77, raising Y p must be compensated for by raising N v in order to avoid 
moving the peak likelihood to higher values of 77 and therefore off of the 7 Li peak. 

In Figure 4, we show the corresponding results with (Li/H) p = 4.1 x 10~ 10 . In this case, 
we must assume that lithium was depleted by a factor of ~ 2.5 or 0.4 dex, corresponding 
to the upper limit derived in ||. The effect of assuming a higher value for the primordial 
abundance of Li/H is that the two peaks in the likelihood function are split apart. Now the 
value of N u = 3 occurs at Y p = 0.227 at i] 10 = 1.1 (a very low value) and at Y p = 0.248 and 
7710 = 5.7. The 95% CL upper limit on N v in this case can even extend up to 6 at Yp = 0.250. 
In Figure 5, we show a compilation of the 95% CL upper limits to N v for different values of 
(Li/H) p = 1.6, 2.0, 2.6, 3.2, and 4.1 x 10 -10 . The upper limit to N„ can be approximated by 
a fit to our results which can be expressed as 

N„ < 80Y P + 2.5 x 10 9 (Li/H) - 15.15 (14) 



6 



Finally we turn to the cases when D /H from quasar absorption systems are also considered 
in the analysis. For the high D/H given in Eq. (f|), though there is formally still a high-?? 
peak, the value of the likelihood function L247 there is so low that it barely falls within the 
95% CL equal likelihood contour (see Figures lb and 2a). Therefore we will ignore it here. In 
Figure 6, we show the peak value of N v and its upper limit for the two cases of (Li/H) p = 1.6 
and 4.1 x 10~ 10 . These results differ only slightly from those shown in Figures 3 and 4. We 
note however, that overall the two values of Li/H do not give an equally good fit. For fixed 
D/H, the high value prefers a value of tjiq ~ 1.8 coinciding with the position of the low-77 
peak for (Li/H) p = 1.6 x 10~ 10 . At higher Li/H, the low-77 peak shifts to lower rj diminishing 
the overlap with D/H. In fact at (Li/H) p > 3.8 x 10~ 10 , the likelihood function L247 takes 
peak values which would lie outside the 95% CL contour of the case (Li/H) p = 1.6 x 1CT 10 . 
The relative values of the likelihood function L247, on the low-77 peak, for the five values of 
Li/H considered are shown in Figure 7. Contrary to our inability to statistically distinguish 
between the two peaks of L47, the large variability in the values of L247 shown in Figure 7 
are statistically relevant. Thus, as claimed in 0], if the high D/H could be confirmed, one 
could set a strong limit on the amount of 7 Li depletion in halo dwarf stars. 

Since the low D/H value has come up somewhat, and since here we are considering the 
possibility for higher values of Y p and (Li/H) p , the statistical treatment of the low D/H case 
is warranted. In Figure 8, we show the peak value and 95% CL upper limit from L247 when 
the low value of D/H is used from Eq. (H) with (Li/H) p = 1.6 x 10~ 10 . The results are not 
significantly different in this case for the other choices of (Li/H) p . In order to obtain N u = 3, 
one needs to go to 4 He abundances as high as Y v = 0.247 with respect to the peak of the 
likelihood function. However, for Y p > 0.234, the revised low value of D/H is compatible 
with 4 He and 7 Li at the 95% CL. The likelihood functions L247 are shown in Figure 9 for 
completeness. 



5 Conclusions 

We have generalized the full two-dimensional (in 77 and N u ) likelihood analysis based on big 
bang nucleosynthesis for a wide range of possible primordial abundances of 4 He and 7 Li. 
Allowing for full freedom in both the baryon-to-photon ratio, 77, and the number of light 
particle degrees of freedom as characterized by the number of light, stable neutrinos, N u , 
we have updated the allowed range in 77 and N u based the higher value of Y v = 0.238 ± 
0.002 ± 0.005 from f2(| which includes the recent data in fi~9fl . The likelihood analysis based 



on 4 He and 7 Li yields the 95% CL upper limits: N u < 4.3 and 7710 < 4.9. The result for N u 
is only slightly altered, N v < 4.4, when the high values of D/H observed in certain quasar 
absorption systems [14", [RJ are included in the analysis. In this case, the upper limit to r/io 



is lowered to 2.4. Since the low values of D/H have been revised upward somewhat [|i~5 
they are now consistent with 4 He and 7 Li and N v = 3 at the 95% CL. We have also shown 
how our results for the upper limit to N v depend on the specific choice for the primordial 
abundance of 4 He and 7 Li. If we assume that the observational determination of 7 Li in halo 
stars is a true indicator the primordial abundance of 7 Li, then the upper limit to N v varies 
from 3.3 - 5.3 for Y p in the range 0.225 - 0.250. If on the other hand, 7 Li is depleted in 
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halo stars by as much as a factor of 2.5, then the upper limit to N v could extend up to 6 at 
Y p = 0.250. 
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Figure 1: (a) L^ 7 (N U , 77) for observed abundances given by eqs. (0 and |3[). 
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Figure 1: (b) L^{N v ,rj) for observed abundances given by eqs. (0, [3], and |). 
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Figure 1: (c) £247(^,77) for observed abundances given by eqs. (|2], |3], and ffy. 
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Figure 2: (a) 50%, 68% & 95% C.L. contours of L 47 and L247 where observed abundances 



are given by eqs. (fj, a, and m. 



13 



5 - 



4 - 



3 - 



2 - 



i — i — r 



J I L 



4- / 

■ : -r ;» < 



t — i — i — r 



D/H = (3.4±0.3)X10 
L 



-5 " 




47 - 
^247 ~Z 



J I I I I L 



10 



•Hio 

Figure 2: (b) 50%, 68% & 95% C.L. contours of L 47 and L247 where observed abundances 
are given by eqs. (0, |3|, and |5|). 
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Figure 3: The position of the value of N v along the low-77 peak (dashed) and high-?] peak 
(dotted) of the likelihood function L 47 as function of Y p . The solid curve shows the 95% CL 
upper limit to N v as a function of Y p . The value of (Li/H) p = 1.6 x 1CT 10 has been fixed. 
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Figure 4: As in Figure § with (Li/H) p = 4.1 x 10" 
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Figure 5: Summary of the upper limits to N v for (Li/H) p = 1.6, 2.0, 2.6, 3.2, and 4.1 x 1CT 10 
as a function of Y p . The lowest curve corresponds to (Li/H) p = 1.6 x 10~ 10 and the limits 
on N v increase with (Li/H) p . 
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Figure 6: As in Figures ^ and ^ based on the likelihood function L247 which includes high 
D/H from quasar absorption systems. 
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Figure 7: Relative values of the likelihood function L247, on the low- 77 peak, for the five 
choices of (Li/H) p in Figure ||. 
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Figure 8: As in Figure |6| for low D/H from quasar absorption systems. 
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Figure 9: As in Figure [7] for low D/H from quasar absorption systems. 
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